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Abstract 

Qh! 

i \ We investigate the corresponding relation between f{R) gravity and an 

interacting holographic dark energy. By obtaining conditions needed for 
c/j ■ some observational evidence such as, positive acceleration expansion of uni- 

verse, crossing the phantom divide line and validity of thermodynamics sec- 
ond law in an interacting HDE model and corresponding it with f{R) mode 
of gravity we find a viable f{R) model which can explain the present uni- 
r~>«. , verse. We also obtain the explicit evolutionary forms of the corresponding 

scalar field, potential and scale factor of universe. 
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1 Introductions 

Observational datafTJ [2j O HI [5], indicates that the current expansion of 
universe is accelerating. Several attempts have been made to justified the 
current accelerated expansion of the universe [6j [11]. One is the 
presentation of an unknown energy form which is called dark energy. An 
alternative approach is the modification of the gravitational theory e.g. f(R) 
gravity in which f{R) is an arbitrary function of the scalar curvature R 
[6j [13l I15| . Recent various observational data imply that the density of 
matter (ordinary matter + dark matter+ radiation), f2 m = 0.27 and the 
density of dark energy, Q\ = 0.73 have capable value today (coincidence 
problem) , beside based on recent data, the equation of state parameter may 
evolve from oo > — 1 (non-phantom phase) in the past to oo < — 1 (phantom 
one) at the present epoch. One way to explain these data, is to consider 
dynamical dark energy with proper interaction with matter [15]. 
In the quantum field theory pa is regarded as zero-point energy density and 
defined based on L, the size of the current universe, (dubbed the holographic 
dark energy ) as follow 

p A = 3c 2 M 2 L" 2 , (1) 

where c 2 is a numerical constant of order unity and M p = l/\/8irG is the 
reduced Planck mass where G is the Newtonian gravitational constant. Dif- 
ferent choices may be adopted for the infrared cutoff of the universe, e.g., 
Hubble horizon, particle horizon, future event horizon [T7]. In a non in- 
teraction model, if we take the particle horizon as the infrared cutoff, the 
accelerated expansion of universe cannot be explained[18j, and if the Hub- 
ble horizon does chooses as the cutoff, then an appropriate equation of state 
parameter for dark matter cannot be derived [19]. By taking, the future 
event horizon as the cutoff, the present expansion of universe may be ex- 
plained but the coincidence problem still remains unsolved. This problem 
may be alleviated by considering suitable interaction between dark matter 
and holigraphic dark energy. 

In this paper, we consider a flat friedmann-Robertson- Walker universe and 
assume that the universe is composed of two interacting perfect fluids, the 
holographic dark energy and the matter. We assume the infrared cutoff 
to be a combination of the future and particle event horizon. After some 
general debate about the properties of the model, we discuss the required 
conditions needed to cross the phantom divide line in the f{R) model. We 
show that this crossing imposes some relations between the parameters of 
the model. 



In this paper we well review f(R) model of gravity and make a corre- 
spondence between f{R) model and an interacting holographic dark energy 
model. By investigating the conditions which are needed for describing the 
present universe, we can obtain a viable f{R) model of gravity. 

2 Description and general properties of the model 

The equation of motion for the f(R) model is 

R^f - \fg^ + (g^a - v^v,) /' = SttGT^, (2) 

where a prime represents the derivative with respect to the curvature scalar 
R and □ is the covariant D'Alembert operator (□ = V a V Q ). We will assume 
dark energy and cold dark matter perfect fluids with stress-energy tensor 
given by 

T iW = -g^vP+ (p + p)u fM u u , (3) 

where p and p are the energy density and pressure of the fluid and u^ = 
(1, 0, 0, 0) is its normalized four-velocity in co-moving coordinates. The dark 
energy component has pressure pd and energy density pa and the cold dark 
matter component has zero pressure and energy density p m . The stress- 
energy tensor is covariantly conserved. 

The trace of equation (J2J) gives an equation of motion for the new scalar 
degree of freedom (compared to Einsteinian general relativity) , [201 El] , 

3D/' = 8ttGT + 2/ - Rf, (4) 

where T is the trace of the stress-energy tensor. It is helpful to redefine the 
scalar degree of freedom through 

<t> = f ~ 1. (5) 

Then Eq. Q can be reexpressed as an equation of motion for a canonical 
dimensionless scalar field <p with a force term T and potential V, 

Ucj) = V'{4>)-F, (6) 

W\4>) = 2f-Rf'l7, (7) 

where the force term that drives the scalar field <j) is proportional to the 
trace of the stress-energy tensor, T = —8ttGT/3. 



Now we consider a homogeneous and spatially-flat spacetime with FLRW 
line element 

ds 2 = dt 2 -a 2 (t){dx 2 + dy 2 + dz 2 ), (8) 

where a(t) is the scale factor. The tt component of the gravitational equa- 
tions d2J, for the metric ([5]), can be simplified to 

The Friedmann equation, @, can be written in a somewhat more conven- 
tional form as 

where we assume that the new scalar degree of freedom behaves like dark 
energy with dark energy density 



3(1 + 
Pd = - 



8irG 



H d(lnf>) lff-Rf 



(11) 



dt 6 V /' 

Also, one can write the Friedmann equation as 

n m + n d = i, (12) 

where the density parameters VL m = p m /p c , &d = Pd/Pci an d the critical 
energy density is 

From the Friedmann equation (|10|) and conservation of stress-energy 
tensor, we have 

H = _^G{ Pt+Pd )_< t H 

1 + 2(1 + 0)' y ' 

where pt = p m + Pd- The vanishing of the covariant divergence of the stress- 
energy tensor for the whole system gives the conservation equation in the 
metric ©, 

pt + 3H( Pt + Pd ) = 0. (15) 

But, because of interactions between the two components, each individual 
component is not necessarily conserved. So, one can write 

Pd + 3H(p d + Pd ) = -Q, (16) 

p m + 3Hp m = Q.in (17) 



We consider different forms of Q below. A number of different models have 
been proposed for dark energy. Here we want to investigate the holographic 
dark energy model and see if it can be related to the f(R) model. Holo- 
graphic dark energy is described in terms of an infrared cut-off length, L, 
and the energy density is defined as 

m Pd = — jpS-. (is) 

where c 2 is a constant of order unity and M p is the Plank mass. This is 
motivated by quantum theory of gravity considerations, in particular the 
holographic principle |18|. fl9l [25] . It was shown in [26] that in quantum field 
theory the UV cutoff A should be related to the IR cutoff L due to a limit 
set by forming a black hole with Schwarzschild radius L. If p d = A 4 is the 
vacuum energy density of the UV cut-off scale, the total energy in volume L 3 
should not exceed the mass of the system-size black hole. This means that 
L s pd < MpL. So for the largest cut-off L, one can define the holographic 
dark energy as (??). From Eqs. (??) and (??), the density parameter of 
holographic dark energy can written as 

The IR cut-off, L, is presumably determined by the available length 
scale. To retain generality, we assume that it is a linear combination of the 
particle horizon, R p , and the future event horizon, Rf, i.e., we choose L to 
be 



where 



L = aRf + fiRp, (20) 



"t-^fwr *' = « ) Lwy <21) 

here t m j n is the time when the particle was created, and ^ a, /3 ^ 1 and 
a+j5 = 1. For a = 1, /3 = we get L = Rf while a = 0, /3 = 1 gives L = R p . 

Taking the time derivative of (??), and using (|16p . one can obtain the 
equation of state parameter uj^ = Pd/pd, 



I22uj d = -1 



1+ Q 209 -a) 



Hp d HL 



(22) 



To progress, we have to specify the interaction term Q. A generic form of 
Q is not available. Three forms which are often discussed in the literature 
are Q = 3b 2 Hpd,3b 2 Hp m ,3b 2 Hpt, where b 2 is a constant which has to be 
positive, because following the second law of thermodynamic, energy transfer 
can only be from dark energy to cold dark mater. These three forms of 
interaction give almost the same result, so for definiteness, we choose 



123Q = 3b 2 Hp d . 
Using OH), (??) and (??), we find 



l2Auj d 



l + 3b 2 



2(/3 



a 



VOTWh 



(23) 



(24) 



As mentioned in the Introduction, observational data indicate that the 
current cosmological expansion is accelerating. In the fluid model this ac- 
celerated expansion requires io d < —1/3. This constraint results in 
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y/(l + <t>)n d < 3b 2 . 



Defining a positive constant < ko < 1, we can rewrite (??) as 

.2(0 -a) 



126- 



V(l + $«o 



3k b 2 . 



(25) 



(26) 



The second law of thermodynamics requires that the entropy S increasing 
with time then, S > 0. We assume S is the entropy attribute to the surface 
area A = A-kL 2 , where L is the infrared cut-off length appearing in (??). 
Also, making use of Nother charge method, one can obtain the entropy in 
the f(R) model of gravity for a horizon with surface A = 4-7T-L 2 as 



Af'(R) 



vrL 2 (l + . 



(27) 



Then, considering the thermodynamics second law, the time derivative of 
entropy, S, should be 
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We set (??) as 
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■SO: 



(28) 



(29) 



where < so- so = is when the accelerating expansion of the horizon of 
universe be adiabatic. Here we assume sq > and then 

, 2(8 - a) 

By making use of (??) and (?? ), we have 

*31 „,/, -, = so - 3fc 6 2 - 2 = O , (31) 

If (1 + 0) 

On the other hand, based on recent data, the dark energy component 
seems to have an equation of state parameter cod < — 1 at the present epoch, 
while ujd > — 1 in the past [23]. Therefore, we expect the equation of state 
parameter cross the phantom divide line, then when cj = — 1, the crossing 
is allowed. So by implying the phantom crossing line constraint on u, (??), 
we have 

Z32 2(/3 ~ a) = 36 2 - 2. (32) 

H L 

From (??) and (??) we have 

Z330 <fco = l-—s<l, (33) 

this show that < 2/36 2 < 1. This means that one of the constant can be 
omit. Moreover, to cross ujd = — 1, &d must be negative at the transition 
time, ujd = — 1. So by using (??) and time derivative of (??) we have 

««,__(* + |tfi) (t»_ |) <a (34) 

From (??), we have b 2 > 2/3, so that the relation (??) is satisfied when 

/35 ;fj > ~\b 2 (35) 

So solving (??), gives 

mH = ^ , (36) 

l + /io7* 

where 7 = 3^o^ 2 /2, Co is an arbitrary constant which satisfy < £0 < 1 
condition. We assume Hoto ~ 1 (Hq and to ar e Hubble parameter in the 
present time respectively) then ho = Hq/(1 — 7). By making use of (??) 



and (??) we can easily find the scale factor of universe and the scalar field 
4> as 

I37a(t) = ao(l + hojt)i , (37) 

<j)(t) = -l + <f> (l + ko'yt)~?,m (38) 

where oo = (1 — 7) ( we assume that the scale factor in the present time 
is equal to 1, a(to) = 1) and 4>q are the integration constants. It is clearly 
seen that, from (??), at the early time a(0) 7^ and then we have bouncing 
in the beginning of the universe. It is well known that the Rcci scalar in flat 
FLRW is as 

J39fl = 6[- + (-) 2 l, (39) 

a a 

then 

where Rq = 6/iq(1 — 27). So that using ©, (??) and (??), one can obtain a 
viable f(R), which allow the crossing from uj = — 1, as 

l41f(R) = f + CR 1 -', (41) 

here Co = (J)qRq 7 /e, e = 6/27 and /o is the constant of integration which 
can be as well as cosmological constant, A. By making use of chameleon 
mechanism this kind of f(R) model has been studied in |27| and they show 
that this form of f(R) model, is viable and satisfy the observational con- 
straints solar system. Also using (??) and (??) one can rewrite (??) as 

IA2V(<P) = V + V 1 4>+ {1 ^ )2ie , (42) 

where Vq is a constant of integration, V\ = 3/o/2 and 

e(l + e)^R 



Vo 



(l-e)(2e-l) 



3 Conclusion 



The HDE model is an attempt for probing the nature of DE within the 
framework of quantum gravity [28]. In this work we used of HDE model 
which is in interaction with DM in the flat FLRW universe. We established 



a correspondence between the interacting HDE model with the f(R) model 
of gravity in the flat FLRW universe. These correspondences are impor- 
tant to understand how different models which have been candidated for 
explanting the present universe, are mutually related to each other. How- 
ever, by taking account an infrared cutoff as a combination of particle and 
future event horizons and using the HDE energy density, we obtained the 
EoS parameter for the interacting HDE. Using equations derived for equa- 
tion of state parameter of holographic dark energy and its time derivative, 
condition required for crossing the phantom divide line was derived. Also 
the condition of validity of thermodynamics second law for the infrared cut- 
off was obtained. Thus we studied the evolving behavior of the interacting 
HDE and implying some observational evidence such as, positive accelera- 
tion expansion of the universe (uj < —1/3 and q < 0), crossing the phantom 
divide line (to < —1) and validity os second law of thermodynamics for an 
interacting model of HDE, we reconstructed the f(R) model which describe 
accelerated expansion of the universe. We also obtained the explicit evolu- 
tionary forms of the corresponding scalar fields, potential and scale factor 
of universe. 
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